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Nuclear  magnetic  resonance  measurements  of  scale  dependent  dynamics  in  a  random  solid  open-cell 
foam  reveal  a  characteristic  length  scale  for  transport  processes  in  this  novel  type  of  porous  medium. 

These  measurements  and  lattice  Boltzmann  simulations  for  a  model  foam  structure  indicate  dynamical 
behavior  analogous  to  lower  porosity  consolidated  granular  porous  media,  despite  extremely  high  porosity 
in  solid  cellular  foams.  Scaling  by  the  measured  characteristic  length  collapses  data  for  different  foam 
structures  as  well  as  consolidated  granular  media.  The  nonequilibrium  statistical  mechanics  theory  of 
preasymptotic  dispersion,  developed  for  hierarchical  porous  media,  is  shown  to  model  the  hydrodynamic 
dispersive  transport  in  a  foam  structure. 


DOI:  1 0. 1 1 03/Phy  sRevLett.  1 03.2 1 800 1 

Solid  open-cell  random  foams  occur  throughout  the 
natural  and  technological  worlds,  ranging  from  natural 
sponges  and  cancellous  bone  to  manufactured  metal,  ce¬ 
ramic,  and  polymer  foams  [1].  In  these  systems,  fluid 
transport  through  the  foam  structure  may  relate  to  biologi¬ 
cal  function  [2]  as  well  as  the  control  of  energy  and  mass 
transport  [3]  in  heat  exchanger  applications  [4,5],  filters, 
catalysts,  and  monolithic  absorbents  [1,6,7].  Character¬ 
ization  of  the  structure  of  foam  systems  is  an  old  and 
important  problem  predating  the  classic  work  of  Kelvin 
and  Plateau  [1,6].  Detailed  geometric  analysis  of  foam 
structure  is  based  on  strut  length,  cell  volume,  and  number 
of  polygonal  faces  and  subsequent  vertex  connections, 
with  ordered  foams  generating  Kelvin’s  tetrakaidecahe- 
dron  and  both  wet  and  solid  random  foams  having  distri¬ 
butions  of  strut  lengths  and  geometric  structures  [1, 6,8,9]. 
Consideration  of  random  foams  as  porous  media  leads  to 
the  concept  of  a  transport  length  scale  [3,10]  and  compari¬ 
son  to  nonequilibrium  statistical  mechanics  transport  mod¬ 
els  [11]. 

In  the  study  of  consolidated  granular  porous  media  the 
length  scale  is  taken  to  be  a  characteristic  pore  size  I  ~ 
(p/{l  —  (p){V/S)  based  on  the  surface  to  volume  ratio 
(5’/y)  and  volume  fraction  cp,  which  derives  from  the 
concept  of  hydraulic  radius  in  fluid  mechanics.  In  highly 
porous  solid  foams,  complications  in  structural  simulation 
and  modeling  are  posed  by  the  solid  phase  forming  a 
sample  spanning  cluster,  which  is  not  the  case  in  lower 
porosity  granular  porous  systems  where  phases  become 
isolated  as  in  percolation  theory  when  P<Pc  [3,12]. 
This  renders  definition  of  a  characteristic  length  scale 
governing  transport  in  foams  an  open  and  important  ques¬ 
tion.  This  Letter  presents  the  first  direct  characterization  of 
the  transport  length  scale  and  nonequilibrium  statistical 
mechanics  modeling  of  transport  in  an  open  cell  solid 
random  foam,  using  nuclear  magnetic  resonance  (NMR) 


PACS  numbers:  81.05.Rm,  05.40. -a,  47.56.+r,  87.64.kj 

measurements  and  lattice  Boltzmann  (LB)  simulation  of 
the  hydrodynamic  dispersion  molecular  dynamics  during 
pressure  driven  fluid  flow. 

NMR  imaging  (MRI)  has  been  used  to  spatially  image 
foam  structures  [13]  and  the  data  used  to  calculate  geo¬ 
metric  properties  [8].  Figure  1  shows  an  optical  image  and 
MRI  of  the  polyurethane  110  ppi  (pores  per  inch,  an 
industry  standard)  open  cell  random  polymer  foam 
(Foamex,  Inc.)  of  volume  fraction  ^  =  0.97  used  in  this 
study.  Alternately,  dynamics  measured  by  pulsed  gradient 
spin  echo  (PGSE)  NMR  provide  length-scale  characteriza¬ 
tion  due  to  pore  structure  sampling  by  molecular  diffusion 
[14,15]  and  advection  [16,17].  PGSE  NMR  measurement 
of  scale  dependent  hydrodynamic  dispersion  has  been  ex¬ 
tensively  applied  to  compacted  granular  media  [17-23]  but 
not  to  characterize  transport  in  open  cell  foams.  NMR 
measurements  [21,23]  and  LB  simulations  [24]  of  hydro- 
dynamic  dispersion  in  random  spherical  bead  packs  dem¬ 
onstrate  a  correlated  motion  in  the  transverse  direction 
perpendicular  to  the  applied  flow.  This  results  in  a  negative 


a)  b) 
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FIG.  1  (color).  Images  of  110  ppi  foam  sample  (a)  optical 
(b)  MRI. 
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FIG.  2  (color).  Velocity  maps  for  the 
longitudinal,  i.e.,  flow,  direction  by 
NMR  for  110  ppi  foam  (sample  1) 

(a)  NMR — transverse  plane, 

(b)  NMR — axial  plane,  and  (c)  LB  simu¬ 
lation  for  50  ppi  foam  of  Montminy  et  al. 
[3] — axial  plane. 


lobe  in  the  transverse  velocity  autocorrelation  function 
(VACF),  and  consequent  preasymptotic  maxima  in  the 
time  dependent  hydrodynamic  dispersion  coefficient, 
analogous  to  that  observed  in  molecular  dynamic  simula¬ 
tions  of  dense  fluids  [25,26]. 

Nonequilibrium  statistical  mechanics  has  been  used  to 
rigorously  derive  a  theory  of  preasymptotic  hydrodynamic 
dispersion  for  porous  media  with  hierarchical  structure 
[1 1]  in  the  context  of  time  correlation  functions  of  dynami¬ 
cal  variables  [25,26].  The  dynamic  variable  of  interest  in 
hydrodynamic  dispersion  is  the  fluctuation  in  Lagrangian 
particle  velocity  u(^)  arising  from  a  Reynold’s  decompo¬ 
sition  of  the  dynamics  into  a  mean  and  fluctuation  compo¬ 
nent  v(^)  =  (v)  +  u(^).  The  evolution  of  the  normalized 
VACF  4f{t)  =  (u(0)  •  u{t)) / {[u(0)y)  is  governed  by  the 
memory  function  equation  [11,25,26] 

=  -  f‘  K(t  -  (1) 

dt  Jo 

A  simple  exponential  memory  function  K{t)  = 
(oloxpi—t/To)  in  Eq.  (1)  generates  the  negative  lobe  of 
the  velocity  autocorrelation  function  [25,26].  The  time 
evolution  of  the  VACF  reflects  the  nonequilibrium  varia¬ 
tion  of  system  transport  due  to  correlations  and  interactions 
of  dynamics.  The  exponential  memory  function  represents 
a  Poisson  process  approximation  to  cage  diffusion  within  a 
fluctuating  cage  in  the  case  of  dense  hard  sphere  fluids 
[25],  and  is  a  model  for  the  position  evolution  of  the 
Brownian  motion  of  a  harmonic  oscillator  [27].  The  expo¬ 
nential  memory  function  is  derived  from  the  second  order 
approximation  of  a  continued  fraction  representation  of  the 
formal  Laplace  transformed  solution  of  Eq.  (1)  [26].  In  the 
hard  sphere  fluid  model,  the  memory  function  parameters 
are  the  second  frequency  moment  of  the  VACE,  i.e.,  the 


force  autocorrelation,  col  =  K(0)  =  (3m)~^V^[/(r)  with  m 
the  fluid  particle  mass  and  [/(r)  the  molecular  energy  po¬ 
tential,  and  T^,  the  correlation  time  of  the  exponential  re¬ 
laxation  to  asymptotic  diffusion  at  equilibrium  [25,26]. 
The  solution  to  Eq.  (1)  for  the  exponential  memory  func¬ 
tion  is 


=  exp(— f/2T^)^cosujr  +  - ^sinujfj, 


(2) 


where  w  =  (co^  —  1/4t^)^/^  and  from  which  the  time 
dependent  hydrodynamic  dispersion  D(t)  =  ([w(0)]^)  X 
/q  is  calculated  from  the  nonequilibrium  Green- 

Kubo  relation. 

Eigure  2(a)  shows  the  spatially  resolved  longitudinal 
velocity  in  the  110  ppi  foam  measured  by  NMR  imaging. 
The  fluid  follows  a  clear  streamline  path  dictated  by  the 
structure  of  the  foam,  with  strong  spatial  correlation  in  the 
flow  direction.  Eigure  2(b)  shows  the  longitudinal  velocity 
in  a  plane  transverse  to  the  flow  direction.  The  velocity 
field  exhibits  little  spatial  correlation  in  the  transverse 
plane.  Eigure  2(c)  shows  longitudinal  velocities  obtained 
by  LB  simulation  of  flow  [28]  through  the  digitized  poly¬ 
urethane  foam  structure  of  Montminy  et  al  [8]  with  50  ppi 
and  ^  =  0.93  (Air  Products,  Inc.).  Although  the  larger  cell 
structure  of  the  50  ppi  foam  is  obvious,  the  similarity 
across  length  scales  between  the  two  foam  structures  in 
Eigs.  2(b)  and  2(c)  is  quite  evident. 

NMR  measurement  of  the  scale  dependent  displacement 
dynamics  transverse  to  the  direction  of  flow  are  shown  in 
Eig.  3  as  a  function  of  the  average  longitudinal  displace¬ 
ment  length  (^||)  =  (u||)A,  obtained  by  varying  the  PGSE 
NMR  displacement  observation  time  A.  Data  were  mea¬ 
sured  for  3  flow  rates  in  the  1 10  ppi  foam  and  based  on  the 
known  linear  scaling  of  transverse  dispersion  with  mean 
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FIG.  3  (color  online).  NMR  measured  hydrodynamic  disper¬ 
sion  coefficient  normalized  by  the  molecular  diffusion  of  the 
fluid  and  the  mean  longitudinal  velocity  (uy)  as  a  function  of 
displacement  observation  time  A  in  terms  of  mean  displacement 
length  (^||)  =  (u||)A.  For  the  110  ppi  foam  sample  1  (open  cir¬ 
cles),  3  flows  ((u||)  =  7.1,  12.1,  and  15.0  mm/s)  were  measured 
and  due  to  linear  scaling  of  with  mean  velocity  have  the 
same  amplitude.  Overlaid  is  the  fit  to  the  diffusion  calculated 
from  the  velocity  autocorrelation  for  the  exponential  memory 
function  model  Eq.  (2)  (blue  or  dark  gray  line).  Data  from  the 
80  ppi  foam  sample  2  at  2  flow  rates  ((uy)  =  6.0  and  8.3  mm/s) 
(open  squares)  and  data  from  the  LB  simulation  (red  or  gray  line) 
for  the  50  ppi  foam  structure  of  Montminy  et  al  [8]  at  (uy)  = 
4.73  mm/s. 

velocity  [21,23,24]  normalization  by  l/(t'||)  yields  a  single 
curve.  All  flow  rates  exhibit  increasing  dispersion  with  a 
maximum  at  the  characteristic  dynamic  transport  length 
scale  I  =  (^||)l(dD^/d<f||»=o  ~  250  induced  by  the 
Structure  of  the  foam  [3,10].  The  data  are  compared  to 
the  hydrodynamic  dispersion  coefficient  calculated  using 
Eq.  (2)  in  the  nonequilibrium  Green-Kubo  relation.  The 
agreement  between  the  theory  and  experiment  is  excellent. 
The  NMR  measurement  directly  provides  the  character¬ 
istic  transport  length  scale  for  the  foam,  thus  complement¬ 
ing  traditional  complex  geometric  analysis  and 
measurements.  Since  the  characterization  and  design  of 
open  cell  foams  is  limited  in  part  by  effective  structural 
characterization  methods,  this  result  is  significant  in  pro¬ 
viding  a  new  and  unique  method  [8]. 

The  transverse  dispersion  for  NMR  experiments  at  2 
flow  rates  in  an  80  ppi  foam  (Eoamex,  Inc.)  also  of  volume 
fraction  ^  =  0.97  and  the  LB  simulation  for  the  50  ppi 
^  =  0.93  foam  of  Montminy  et  al  [8]  are  shown  in  Eig.  3. 
The  dispersion  behavior  is  qualitatively  similar  in  all  three 
cases  but  the  exact  location  of  the  maximum,  the  transport 
length,  depends  on  the  porous  structure. 

Eigure  4(a)  demonstrates  data  collapse  with  the  dis¬ 
placement  length  (^11 )  scaled  by  the  transport  length  I  and 
the  dispersion  coefficient  normalized  by  the  maximum  in 
Dx((^l|))  for  the  foam  systems.  The  transport  lengths  are 
I  ~  250  jmm  for  the  1 10  ppi  foam  (sample  1),  I  ~  300  /mm 


for  the  80  ppi  foam  (sample  2)  and  /  ~  475  /mm  for  the 
50  ppi  foam  [8].  The  nonequilibrium  statistical  mechanics 
model  scales  identically,  in  excellent  agreement  with  the 
NMR  data  and  LB  simulation.  Also  plotted  is  LB  simula¬ 
tion  of  a  random  sphere  pack  of  ^  =  0.44  [24]  in  which  the 
transport  scaling  length  used  is  /  ~  ^/(l  —  cp)d  =  O.^d 
for  spheres  of  diameter  d.  Remarkably,  despite  significant 
structural  and  porosity  differences,  the  sphere  data  collap¬ 
ses  with  the  foam  data  for  I  <  (^y),  and  shows  only  slight 
variation  at  longer  displacements.  The  consistency  be¬ 
tween  these  results  and  NMR  data  and  LB  simulations 
for  consolidated  spherical  grain  porous  media  of  much 
smaller  porosity  [21,23,24],  indicates  the  universality  of 
dynamics  over  a  broader  range  of  porous  structures  than 
previously  recognized.  Converting  the  hydrodynamic  dis¬ 
persion  into  a  variance  and  taking  the  second  time  deriva¬ 
tive  of  the  NMR  data  for  each  flow  rate  of  the  110  ppi 
foam,  then  averaging  the  amplitude  scaled  values  yields 
the  velocity  autocorrelation  data  in  Eig.  4(b).  Again,  when 
plotted  versus  longitudinal  displacement  scaled  by  the 
transport  length,  the  agreement  between  experiment,  simu¬ 
lation  and  analytical  theory  is  excellent. 

The  results  of  Eigs.  3  and  4  indicate  the  potential  for 
further  development  of  memory  function  approaches  to 
model  transport  in  porous  media.  The  ability  to  incorporate 
pore  structure  into  model  parameters  for  the  scale  depen¬ 
dence  of  dynamics  provides  a  means  to  design  porous 
structures  which  control  transport.  In  the  exponential  mem¬ 
ory  function  model  for  dispersion,  the  col  represents  a 
potential  which  determines  the  transverse  velocity  oscilla¬ 
tions  due  to  transport  around  the  foam  structure,  analogous 
to  the  strength  of  the  restoring  force  in  the  harmonic 
oscillator  analogy  [27]  and  consistent  with  an  oscillatory 
flow  toy  model  employed  to  elucidate  the  physics  of  simi¬ 
lar  dynamics  in  consolidated  granular  porous  media  [21]. 
Eollowing  the  arguments  of  Berne  et  al,  the  negative  lobe 
of  the  VACE  in  Eig.  4(b)  indicates  reversal  of  fluid  motion 
due  to  the  foam  structure  [25].  The  mean  displacement  a 
fluid  particle  with  initial  velocity  u{0)  travels  before  ve¬ 
locity  reversal  (AY(^))  =  jl{u{t'))dt'  =  J  u{0)il/{t')dt'  is 
given  by  the  first  zero  crossing  of  the  VACE  [25].  In  a  pure 
oscillatory  flow  such  as  the  toy  model  of  Callaghan  and 
Codd  [21]  this  length  is  the  order  of  a  quarter  wavelength 
of  the  oscillation.  The  transport  length  scale  I  defined  by 
the  maxima  in  the  transverse  dispersion  corresponds  to  the 
length  scale  at  the  minima  of  the  VACE  and  is  half  the 
wavelength  of  a  pure  oscillatory  flow.  In  the  foam,  this 
transport  length  represents  the  mean  intercell  distance  over 
which  a  fluid  particle  experiences  the  entire  transverse 
velocity  distribution.  The  heterogeneity  of  the  velocity 
field  induces  a  decay  of  this  correlation  with  time  scale 
T^.  Development  of  a  potential  function  based  on  the 
porous  structure  that  generates  particular  fluid  path  lines 
would  represent  a  major  advance  toward  connecting  struc¬ 
ture  and  transport  in  porous  systems  [3,12]. 

Hydrodynamic  dispersion  of  fluid  flow  through  open  cell 
foam  has  been  demonstrated  to  provide  a  dynamic  trans- 


218001-3 


PRL  103,  218001  (2009) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
20  NOVEMBER  2009 


FIG.  4  (color  online),  (a)  The  transverse  dispersion  coefficient  normalized  by  its  maximum  amplitude  as  a  function  of  the  axial 
displacement  length  scaled  by  the  transport  length  scale  /  for  LB  simulation  (red  or  gray  line)  of  the  foam  structure  of  Montminy  et  al. 
[8],  the  NMR  data  for  the  1 10  ppi  foam  sample  1  (open  circles)  at  (vy)  =  12. 1,  the  80  ppi  sample  2  (open  squares)  at  (vy)  =  8.3  mm/ s, 
the  scaled  memory  function  model  (blue  or  dark  gray  line),  and  LB  simulation  for  a  consolidated  random  sphere  packing  24  (black 
line),  (b)  The  velocity  autocorrelation  for  the  average  of  all  the  NMR  experiments  on  the  110  ppi  foam  sample  1  and  the 
nonequilibrium  statistical  mechanics  model  (blue  or  dark  gray  line)  and  L-B  simulation  (red  or  gray  line). 


port  length  scale  based  on  the  correlated  dynamics  of  the 
fluid  transverse  to  the  flow  direction.  This  directly  provides 
a  characteristic  length  scale  for  the  foam  without  geometric 
analysis.  The  transverse  hydrodynamic  dispersion  dynam¬ 
ics  measured  by  NMR  for  two  different  foam  samples,  LB 
simulation  for  a  third  foam  and  random  sphere  pack,  and 
the  analytic  memory  function  model  are  shown  to  collapse 
when  scaled  by  this  transport  length  scale  and  the  maxi¬ 
mum  amplitude  of  the  transverse  hydrodynamic  disper¬ 
sion,  indicating  a  universality  of  the  dynamics. 
Agreement  between  nonequilibrium  theory,  simulation 
and  experiment  is  established  providing  a  basis  for  further 
study  of  a  broad  range  of  porous  materials. 
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